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Abstract — The study of the scattering data for a star-shape 
network of LC-transmission lines is transformed into the 
scattering analysis of a Schrodinger operator on the same 
graph. The boundary conditions coming from the Kirchhoff 
rules ensure the existence of a unique self -adjoint extension of 
the mentioned Schrodinger operator. While the graph consists 
of a number of infinite branches and a number finite ones, 
all joining at a central node, we provide a construction of 
the scattering solutions. Under non-degenerate circumstances 
(different wave travelling times for finite branches), we show 
that the study of the reflection coefficient in the high-frequency 
regime must provide us with the number of the infinite 
branches as well as the the wave travelling times for finite 
ones. 

I. Introduction 

The number of electronic equipments is increasing 
rapidly in automotive vehicles, aircrafts, and many other 
safety critical systems. Consequently, the reliability of wired 
networks and electric connections is becoming more and 
more important. For example, in automotive industry, a 
goal is to develop compact and easy to use devices for 
the diagnosis of electric connection failures in garage or at 
the end of the production chain. These devices should be 
capable of detecting and locating failures in cables and in 
connectors. Another goal is on-board diagnosis: the diagno- 
sis device will be integrated to the vehicle in order to detect 
failures under normal working conditions of the vehicle. 
To find faulty wiring in such networks, it is not always 
possible to measure end-to-end cable impedances, because 
the number of available diagnostic port plugs is limited, 
and furthermore, for diagnosis purpose, it is not sufficient 
to detect a high end-to-end impedance as it is also necessary 
to locate the fault within the cable. In such situations Time 
or Frequency Domain Reflectometry (TDR, FDR) are the 
most commonly used methods: a high frequency signal (a 
short travelling pulse for TDR, a standing wave for FDR) is 
sent down a wire at some point and the signal reflected by 
the network is measured at the same point and analyzed for 
fault detection and location [11], [16], [5], [2]. Automatic 
fault detection and diagnosis using reflectometry methods 
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is the subject of intense research, both on the technologies 
of "smart wiring systems" and reflectometers [15] and on 
the foundations of the TDR/FDR methods. Most studies are 
concerned with various aspects of mathematical modeling 
and simulation of microwave propagation in networks, from 
network complexity [12] to reduced order modeling of 
high frequency phenomena like skin effect [1]. The inverse 
problem of fault detection and localization for a network 
is in general studied in simulation or experimentally. For 
this problem, the available theory is developed mainly in 
the framework of Inverse Scattering Theory limited to very 
simple networks (segment, half-line, line) (see e.g. [3], 
[13]). 

In this paper, we consider the inverse scattering problem 
(ISP) for a star-shape network F of transmission lines. Some 
of the branches of the network are extended in a direction z 
to +°° and the others admit a finite length. The network is 
assumed to be non-dissipative; i.e. the series resistance R(z) 
and the shunt conductance G(z) per unit length are assumed 
to vanish and therefore, we are dealing with LC-wires. 

The main concern here is to derive some information 
concerning the topology and the geometry of the network, 
applying a very few scattering information. Indeed, fixing 
an infinite branch e\ (see Figure [TJ, we generate a high- 
frequency time varying voltage at its infinite end and 
measuring the intensity at the same end, we would like to 
derive information such as the number of infinite branches 
m and the wave travelling time z m+ j of the finite segments 

Cm+J- 

Here, in order to model our network, each branch is 
parameterized through its own length. The central node is 
assumed to be the zero (origin) of all the branches. While 
the position coordinate Zj takes values inside the interval 
[0,lj], lj being the length of the branch <? ; - and therefore 

lj = °° Vje{l,---,m} and lj<°° V/ € {m+ 1, • • • ,m+n}. 

Moreover, on each branch, of finite or infinite length, we 
suppose that 

Al L(z) the inductance and C(z) the capacitance are 

sufficiently regular (twice differentiable); 
A2 L(z) > and C(z) > 0. 

A3 On infinite branches, L(z) and C(z) have strictly 
positive finite limits and C(°°) as z — > °°. 

In the next section, we provide the mathematical model 
consisting of the transmission line equations on each one of 
the network branches and the boundary conditions (Kirch- 
hoff rules) coupling these equations together. Next, for the 




Fig. 1 . The LC-network T consisting of m infinite branches extended to 
Z = °° and n finite branches of various lengthes. 



case of a simple transmission line, we provide a brief 
review of an old result by I. Kay [9], transforming the 
inverse scattering problem for one-dimensional nonuniform 
non dissipative transmission lines to the inverse scattering 
problem for a Schrodinger operator. We extend this result 
to the case of our network by defining relevant boundary 
conditions. Finally, we provide a result ensuring the es- 
sential self-adjointness of the Schrodinger operator on the 
graph with the aforementioned boundary conditions. The 
existence of a unique self-adjoint extension allows us to 
start the study of the scattering theory and the scattering 
solutions. 

Next, in Section [Till we solve the scattering problem by 
providing the scattering solutions. In this aim, we apply 
the Jost solutions on semi-infinite branches [7], [6] and the 
fundamental solutions on the finite segments [10]. 

Finally, in section |IV] we study the asymptotic behavior 
of the reflection coefficient in the limit of high-frequency 
waves and we prove how this limit provides us the needed 
information on the topology and the geometry of the graph. 

II. Mathematical model 

A. Transmission line equations 

In the absence of dissipation, the transmission line equa- 
tions writes 



di 



(i) 



where I(z,t) and U(z,t), respectively denote the intensity 
of the current and the voltage at time t and position z of 
a finite or infinite branch. Let us index these functions as 
{Lj}"jt"\ {CjYjti, {IjYjti and {Uj}^, defined on the 
branches {e;}"±f . 

At the central node, the Kirchhoff rules imply 



Ut(0)=Uj(0), 

m-\-n 

£/;(0)=0. 



ij 6 {!,-•• ,m + n} 



At the terminal nodes l m +j of the finite branch e 
assume the boundary condition 

l m+ j(l m+ j) = ;£{!)•" )«}• 



(2) 



m+j, we 



(3) 



B. Liouville transformation 

Here, we follow an old result by I. Kay [9] and extended 
by M. Jaulent to the dissipative case [8]. 
For a wave of frequency k, i.e. for 

I{z,t)=I(k,z)e-' kt 
U(z,t) = U(k,z)e- ,kr , 

the equation (fTJ becomes 

d ( 1 dU" 



dz \L(z) dz 
Applying the Liouville transforation 



k z C(z)U = 



(4) 



(L(m)C(m)) 2 du 



and the conventions I(k, z(x)) =I(k,x), U(k,z(x)) = U(k,x) 
etc, and setting 



y(k,x) 

the equation writes 



C(x) 
L(x) 



U(k,x) 



(5) 



^y + (k 2 -V(x))y = 0. 
Here, the potential V(x) is defined as follows: 

V(x) 



\C(xY 


- * d 2 


[C(x)l 




dx 2 


[L(x)\ 



(6) 



Note that, the new coordinate x here has the dimension 
of time. It actually parameterizes the travelling time of 
the wave through the transmission line. In particular the 
length / of a finite transmission line is transformed into the 
correspoindg wave travelling time: 



= f {L{z)C{z)) l l 2 dz 
Jo 



C. Extension to the network 

The LC-transmission line equations defined on each edge 
for j e {l,...,n +m} are equivalent to the Schrodinger 
equations 

d 2 yj 



dx 2 



+ (k 2 -V j {x))yj = 0. 



(7) 



In this new formalism, the boundary condition ([2j writes 

■y ij G {I,-" ,m+n}, 

(8) 



A 7V/ 



v=0 



x=0 



Here, applying the assumptions of the Section |H the func- 
tions 



Aj( X j) 



Cj(xj) 



j G {!,■■■ ,n + m} 



- L A x .i). 

are well defined, strictly positive and regular. 



Furthermore, the boundary conditions (0 write 



III. Scattering Solutions 



ym+ji^m+j) ■ 



A' 



m+j 



1 m+ j 1 



j 6{l,-,n}, (9) 



where {%n+jYt=\ denote the wave travelling time corre- 
sponding the finite transmission lines. 
In this paper, we denote 



A' 



1m+j 



m+j 



*m+ j 



j =!,■■■ ,n. 



(10) 



In conclusion, in order to study the LC-transmission line 
equations on the graph T with boundary conditions (O 
and (01, we can study the Schrodinger equation dT) with 
boundary conditions (JHJ and (O. 



We are interested in the scattering solution where a signal 
of frequency k is applied at one end of one of the infinite 
branches. In fact, in general, one might be interested in a 
experiment where signals of different amplitudes {dCyV-Lj 
are applied at all ends of the infinite branches. In such a 
case, we will be seeking a scattering solution satisfying the 
asymptotic behavior at k — > °° 



y J {x j ,k)~a j e lkx > 



-Rj(k)e l 



k.Xj 



However, the experiment we consider here corresponds to 
the simple case 

a\ — 1 and a,- = 0, j = 2,--- ,m. 

The main result of this section, may be expressed as follows: 



D. Self-adjointness 

We are dealing here with a symmetric operator Jz? defined 
onL 2 (r): 



dxj 



J/? = ®5?j, £?j = -^-2+V(xj) on ej, ;'=!,••• ,m+n, 



with 



D(j£?) = closure of C£ c in H 2 {T). 



Here, C£° c denotes the space of infinitely differentiable 
functions satisfying the boundary conditions <[8J an d ®- 

In order to start the study of the scattering theory for 
this operator, we need first to prove that it admits a self- 
adjoint extension. In this aim, we assume an appropriate 
decay assumption on V(xj) ensuring that the operator Jz? 

is a compact perturbation of the operator ^— J^- 

Indeed, one also need to add a new assumption A4 to the 
three assumption Al through A3 of the introduction: 

A4 On each branch ej, j — 1 , ■ ■ ■ ,m + n, we have 



(1 



c A x j) 
L j( x j) 



72 r 



dx) 



<^ oo. 



Now, we apply a general result by Carlson [4] on the self- 
adjointness of differential operators on graphs. 

Indeed, following the Theorem 3.4 of [4], we only need to 
show that at a node connecting /V edges, we have /V linearly 
independent linear boundary conditions. At the terminal 
edges of this is trivially the case as there is one 

branch and one boundary condition. At the central edge it 
is not hard to verify that ([H) and (0 define such boundary 
conditions as well. This implies that the operator Jz? is 
essentially self-adjoint and therefore that it admits a unique 
self-adjoint extension. We are now ready to start the study 
of the scattering solutions. 



Theorem 1: There exists a unique solution 

yM) = (y,(*^));t + ;, 

of the scattering problem, satisfying 

. if/y j (xj , k) = k 2 y j(xj,k); 

' (yj( x jik))'j=" satisfy the boundary conditions ^ 
and ©; 

• For each k £ R, there exist Ri(k) and {Tj{k)} k : =2 such 



that 



,/A.vi 



yj(xj,k)~Tj(k)e lkx 3, a i 



xi^°°, (11) 
j = 2,---,m. (12) 



The coefficients, R\ (k) and {Tj}J =2 , in (fTTT i and (TTZb appear 
to be unique. They are called respectively the reflection and 
the transmission coefficients. 

We proceed the proof of the Theorem Q] by constructing 
such a solution and by showing that such a construction is 
unique. 

In this aim, we apply two different types of solutions: 
1- the so-called Jost solutions on half-lines [7], [6] and 2- 
the fundamental solution of the one-end Sturm-Liouville 
equation [10]. 

A. Jost solutions 

The Jost solutions f(x,k;V) and f(x,k;V) on the half 
line [0,°°), where V is a potential defined on [0,°°), are the 
solutions of the Schrodinger integral equations 



f(x,k;V) 
f(x,k;V) 



sink(x — y) 



V(y)f(y,k;V)dy, 



-ikx 



sink(x — y) 



V(y)f(y,k;V)dy. (13) 



Under the decay assumption A4 on V, these Volterra type 
integral equations admit unique solutions. In fact, one only 
needs to apply the standard iterative methods for the proof 
of existence and uniqueness of solutions to Volterra type 
integral equations (see e.g. [14], pages 138-139, for a proof). 



The functions f(x,k;V) and f(x,k;V) satisfy the the 
Schrodinger equation 

-& f + V{x)f = k2f > -&2f + V ( X )f = ^f 

and admit the asymptotic behaviors (at x — > °°) 

f(x,k;V)=e lkx + o(l), 
f(x,k;V) =a{k;V)e-' kx + b{k;V)e ,kx + o{l) 7 

where 

a(k;V) = 1 -^J°°e lkx V(x)f(x,k;V)dx, 



and 



(14) 



b(k;V) 



1 



- (fee 



cx V(x)f(x,k;V)dx. 



Iik Jo 

Moreover, one can show that ([14], pages 138-139) 
\^f(x,k;V)-ike lkx \<\k\\^p(^f\V(y)\dy)-l\. (15) 

B. Fundamental solution 

The fundamental solution, co(x,k;h,V), is a solution of 
the Schrodinger equation, 



a j 

--rffllVWffl = 0) 



ra(0,ifc;fc,y) 



1 



*e(o,»)> 

fl)'(0 1 jfc;A i y)=A, 



(16) 



where ©' denotes the derivative with respect to x. 

This solution may be expressed through the following 
integral representation [10], 



<x)(x,k;h,V) = cos(fac) +h 



sin(fac) 



+ / K(x,t;V){cos(kt)+h 



sin (kt ) 



}dt. (17) 



Here the kernel K(x,t;V) is a solution of the following 
integral equation 



K(x,t;V) 



1 

2 Jo 



V(u)du 



l f T~ 



rfa / y(a + j3)^(a + j3,a-j3)Jj3, (18) 
2 Jo Jo 

where we have extended the potential V to be zero outside 
the interval of its definition. 

The integral equation ([T8l admits a unique solution 
(cf. [10], Theorem 1.2.2). Applying the same result, for a 
potential V in L (M), this solution satisfies 

\K(x,t;V)\ < — IIVH^i expCjcHVHii), for x > and \t\ <x. 



Furthermore, simple computations imply 

K x (x,t;V) = -V( — ) 



(19) 



4 Jo 



[ 2 da V(a- 

10 



x-t 



x — t x — t. 
-—^-— ] 



1 ,x + t, 
K t (x,t;V) = -V(—) 

1 



- x+i 

/ daV(a+ )K(a- 

JO 2 



4,/o 



X — t x-t, 



where K x and K t denote respectively the derivatives with 
respect to the first and the second coordinates x and t. For 
a potential V G L 1 (R), this together with $1% implies 

\K x (x,t;V)\,\K t (x,t;V)\<^\V(^)\ 

+ U\V\\ 2 Ll exp(x\\V\\ L i), forx>0 and \t\ <x. (20) 

C. Proof of Theorem\l} constructing the scattering solution 

For this construction, we distinguish between three types 
of edges. 

1) Edge e\: Define 

1 



C(x u k) =yi(xi,k) 



a(k;Vi) 



Vi). 



Note that 



a(k;V) = l + — / V(x)dx + o(l/\k\), (21) 
2ik Jo 

and therefore, at least for k large enough, a(k,V\) is non- 
zero and £(x\,k) is well-defined. Furthermore, note that 
£(xi,k) satisfies the Schrodinger equation 

-^;+v 1 (x 1 ); = k 2 ;, 



and through the requirement (TTTT) and by (TPfl i. at x\ — > 00 

(22) 



a(fe;y) 



Consider now the Wronskian 



W(a,k)J(.,k;Vi))(xi) = 

C (jc, , (jti , k; Vi) - C'Oci , Vi) , 

where the derivatives are taken with respect to the position 
X\. It is not a hard task to verify that the above Wronskian 
is a constant function of x\. In fact, as the both functions 
and f(.,k;Vi) verify the same Schrodinger equation 
with the same potential V\, one can easily check that 

^-W(C(.,k),f(.,k-V l ))(x l )=0. 
dx\ 

Applying ( l22l and ( TT3T >. we have 

lim W(Z(.,k)J{.,k;Vi))( Xl )=0 

and therefore 

W(£(.,*),/(.,*;Vi))(jci) =0 for ^ e [0,00). 



This implies that £(.,&) and f(.,k;V\) are co-linear and thus 
for some R\ (k) 



yi(xi,k) 



1 



a(k;Vi) 



+ (Ri(k)-^\)f(xx,k;Vx). (23) 

2) Edges ej, j =2, ■■■ ,m: We consider the Wronskian 
W(yj(.,k),f(.,k;Vj))(xj). Just as in above, this Wronskian 
is a constant of xj and applying (fT2l together with ([T5"l >. 

lim W(yj(.,k)J(.,k;Vj))(xj)=0 



and thus for some Tj (k) 

yj(xj,k) = Tj(k)f(xj,k;Vj), 



]=!,■■■ ,m. 



(24) 



3) Edges ej, j = m + !,■■■ ,m + n: We consider the 
Wronskian 

W(jj(Tj - .,k), co(tj - .,k;hj,Vj(Tj - .)))(*;), 

where hj is defined in (TTOb . On the segment [0,T/], the 
functions yj(?j — -,k) and (o{ij — .,k;hj,Vj(Tj ~ .)) satisfy 
the same Schrodinger equation with the same potential 
Vj(tj — .), and therefore the Wronskian is a constant. 

Applying (0 together with the definition of the funda- 
mental solution ( D~6l >, we have 



W(yj{.,k),tofa-;Mj,Vj(Tj-.))) 



= 0, 



and therefore, just as in above there exists a constant ttj(k) 
such that 

yj(xj,k) = aj(k)co(Tj -Xj,k;hj,Vj(Tj - .)), 

j = m+ l,--- ,m + n. (25) 

At this point, we can apply the boundary conditions at the 
centeral node ([§), in order to find the m + n constants R\(k), 
{Tj{k)}J =2 and {cCj(k)}J+^ +v In fact, we have 



y=V(o) 



1 



a(k;Vi) 



(*i(*)-^r^)/(o,*;Vi) 



a(k;Vi) 



-A^(0)T h {k)f(0MV h ) 



(26) 



for j\ = 2, • • • ,m and j'2 = m+ 1, ■ • • , m + n. Here, for the 
first equality, we have applied f(0,k\V\) — 1. Moreover, we 
have 

a(k\V\) a[k;Vi) 

in 

+ ^Aj(0)Tj(k)f'(0,k;Vj) 

7=2 

m+« m+n 

+ I AjiftaiWal^kihMxj - .))=y £ A;(0)4(0). 

;'=m+l 7=1 

(27) 

Here, for the first line we have applied f'(Q,k;V) = ik. 
One can easily check that, ( |26| i and ( |27] i provide m + n 



linearly independent linear equations for the m + n unknown 
constants Ri(k), {7}(A;)}™ =2 and {<*j}1 =m+ \ - One can there- 
fore find these coefficients uniquely and solve the scattering 
problem. □ 

Here, we are only interested in the high frequency (k — » 
°°) behavior of the reflection coefficient R i (k) and we do 
not need to compute all the coefficients explicitly. 

IV. High-frequency asymptotic behavior of 

REFLECTION COEFFICIENT 

In this Section, we announce and prove the main result 
of this paper: 

Theorem 2: Consider the assumptions Al through A4, 
together with the new one, 

A5 The function L/C is continous at the central node 
od T; i.e. 



LjM _L h (0) 



,m + n. 



c A (o) c h (oy 

Then, in the high-frequency regime (k — » °°) the reflection 
coefficient R i (k) satisfy: 



Rl(k) 



(m-2>-^ ±1 tan(^ 
™-I-L + „" + itan(^) 



+ 0(1). (28) 



Note that, applying this result one can easily identify m 
the number of the infinite branches, as well as {T ; }™j JJ " +1 
the wave travelling times of the finite ones (when they 
are all different) through the reflection coefficient R\(k). 
Indeed, through simple computations, we have the main 
result of this paper: 

Corollary 1: Under the assumptions Al through A5, and 
in the high frequency regime k — > °°, we have 

m = 29t( - ■ I ,„ )+o(l), 



(29) 



and 



(30) 



7/7 -h 77 1 

£ tan(H / ) = 23(— — -)+ (l) 

j=m+l l+Ri(k) 

where 3i and 3 denote respectively the real and imaginary 
parts of a complex number. In particular, scanning a certain 
interval of frequencies k for k > k* large enough, and 
through the study of the poles of 3( l+ l R ^ ), we can identify 
the wave travelling times Tj. 

Before starting the proof of the Theorem [2] we need to 
prove two lemmas. 

Lemma 1: Under the assumption A4 on V(x), we have 

f'(0,k;V) 



f(0,k;V) 



ik+0{\) ask 



(31) 



Proof of Lemma [7} Here, we just need to apply the 
the following approximations on f(x,k;V) and f (x,k;V) 
(cf. [14], pages 138-139, for a proof): 



\f(x,b,V)-e' kx \<\eM±[\V(s)\ds) 
,df(x,k;V) 



- 1 



dx 



-ike lKX \ < \k\\exp( 



\k\ 



V(s)\ds)-l\ 



(32) 



Applying the integrability of V(x), the relation ( BTT i can be 
obtained very easily. □ 

Lemma 2: Under the assumption A4 on V(x) and for a 
fixed position x = T, we have 

m'(r k-h V\ 

) 1 ,1, - ~-fctan(feT)+o(l) asfe^oo. (33) 

Proof of Lemma^ We apply here the expression dTvT > of 
0)(x,A:;/z, V). Through the inequality dT9l >, we now that for 
a fixed T, the kernel ^T(t,?;V) is bounded and therefore we 
easily have 



reflection coefficient may be derived through the Lemmas Q] 
and |2] and the relations (T2T1 i and 



(0{z,k;h,V) =cos(£t) + J K(T,t;V)cos(kt)dt + 0(l/k). 
developing the integral by parts, we have 

/ K(T,t;V)cos(kt)dt= ~sin(to)K(T,t;V) ' 

J-T K t=-T 

1 f T 

/ sm(kt)K t (x,V,V)dt 

k J-t 

Now, applying the inequality d20l >. we have 
/ K(z,t;V)cos(kt)dt < -sin(kt)K(z,t;V) 

J-T k t=-T 

where we have applied the boundedness of K(x,t;V) for 
\t\ < T and the integrability of V. This implies 

(0(T,kr,h ) V)=cos(kT) + 0(l/k), as£^°°. (34) 
Furthermore, we have 

a x {T,k;h,V) = -ksin(kT) + hcos(kx) 



+ (K(r, t) + K(%, -t)) cos(fcT) + /j (A:(t, t) - £(t, -t)) 



sin(fe) 



+ K x (z,t;V){cos{kt) + 



sin(fa) 



}dt. 



Similar computations as in above, together with the inequal- 
ity (1201 . imply 

(0 x (T,k;h,V) = -ksm{kT)+O(l), ask^™. (35) 

The two relations (134-b and (f35T > finish the proof of Lemma [2] 
and we have (l33l l. □ 

We are now ready to prove the Theorem [2] 

Proof of Theorem [2} 

The assumption A5 implies the existence of a constant 
A, such that 

A 7 -(0)=A, \/j =!,■■■ ,m + n. 
Dividing (|27T > by y and applying (l26l i, we have 

"Ya' (0) I A = — + {a{t ' Vl )Rl {k) - ^ Vl ] )f ( °' k Vl ] 

p x 1 i + (a(*;Vi)iii(*)-*(fc:Vi))/(o,k;Vi) 

- /(0,fe;Vy) m +* o/^.fci^.VQ-^-.)) 

je 2 /(o,*:^) J= tVi co(Tj,k;hj,Vj(zj - .)) ■ 

Through this relation we can derive the value of the reflec- 
tion coefficient Ri(k). The high-frequency behavior of this 



b(k;V) < \\V\\ L iO(l/k) ask 



(36) 



Some simple computations, then imply the relation (l28l i and 
finish the proof of the Theorem |2] □ 

V. Conclusion 

As we have seen, the high-frequency asymptotic behavior 
of the reflection coefficient, over a fixed semi-infinite edge 
of a star-shape graph, must provide us very interesting and 
useful information on the geometry and the topology of 
our network. Two main directions may be considered for 
further extension of this result. The first direction deals with 
the more general case of trees instead of simply star-shape 
graphs. The second one corresponds to the case of LCRG 
transmission lines where one needs to extend all the utilities 
to the Zakharov-Shabat operator [8]. 
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